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Abstract 

Based on the properties of probability distributions of functions of random vari- 
ables, we proposed earlier a simple stringy mechanism that prefers the meta-stable 
vacua with a small cosmological constant A. As an illustration of this approach, we 
study in this paper particularly simple but non-trivial models of the Kahler uplift in 
the large volume flux compactification scenario in Type IIB string theory, where all 
parameters introduced in the model are treated either as fixed constants motivated 
by physics, or as random variables with some given probability distributions. We 
determine the value Wq of the superpotential Wq at the supersymmetric minima, 
and find that the resulting probability distribution P(wq) peaks at wo = 0; further- 
more, this peaking behavior strengthens as the number of complex structure moduli 
increases. The resulting probability distribution -P(A) for meta-stable vacua also 
peaks as A —7- 0, for both positive and negative A. This peaking/ divergent behavior 
of -P(A) strengthens as the number of moduli increases. In some scenarios for A > 0, 
its expectation value (A) decreases as the number of moduli increases. The light 
cosmological moduli problem accompanying a very small A is also discussed. 
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1 Introduction 

Recent cosmological data strongly suggests that our universe has an exponentially small positive 
cosmological constant A ~ 10~ 122 Mp (see [I] and references therein). On the other hand, string 
theory has so many possible meta-stable vacuum solutions that it should have at least a solution 
with such a small A |2|. So it leaves open the question why nature picks a vacuum solution with 
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such a very small A, in units of the Planck scale. In a previous paper [3], we propose a plausible 
reason why this may happen within the context of string theory. 

The basic idea is very simple. A typical compactification in string theory involves many 
moduli and fluxes (see the review [4]). The moduli and their dynamics describe the string 
theory landscape. Stabilizing them will lead to a set of values for these moduli, so A of a meta- 
stable vacuum will be a function of them. Many if not most of the moduli can take multiple 
(some discrete) values since each stabilized compactification involves a set of quantized fluxes. 
As a result, many of these moduli will take values within a range that include zero. If (and this 
is a big if) A is a product of some of them, then the probability distribution -P(A) of A will 
naturally peak (even mildly diverge) at A = 0. In actual models, the functional dependence of 
A on the parameters in the model is much more complicated. Here we show that the peaking 
property of -P(A) is present in the stringy scenarios studied here. 

In [3], we work out the case for a single Kahler modulus in the Kahler uplifting region |5}j8] of 
the Large Volume Scenario |9] in Type IIB string theory. There we treat the few parameters in 
the model as random variables and show that a modest suppression of A is achieved: despite of 
the non-trivial functional dependence of A on the parameters, -P(A) is peaked (actually diverges) 
at A = 0, even though the parameters are treated as random variables with uniform (or similarly 
smooth) distributions. In this single Kahler modulus model [3], we see that the suppression of 
A is present but very modest. In this paper, we consider the multi-moduli cases to check the 
validity of the basic idea. Here our probability and statistical analysis is largely based on the 
study by Rummel and Westphal [7] on a simplified yet non-trivial model in the Large Volume 
Scenario. Although the functional dependence of A on the parameters is more non-trivial, we 
still find that the probability distribution -P(A) for meta-stable vacua becomes more peaked 
at A = (even mildly diverges) as the number of moduli increases. However, whether the 
expected value of A decreases exponentially or not as the number of moduli increases depends 
on the details of the physics which remains to be better understood. 

To be more specific, consider a Swiss-cheese type of Calabi-Yau three-fold with h 1 ' 1 number 
of Kahler moduli and h 2,1 number of complex structure moduli (so the manifold M has Euler 
number x(M) = 2(/i 1 ' 1 - h 2 ' 1 ) and we are mainly interested in negative x(M)). The simplified 
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model of interest is given by, setting Mp 
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where Q is the homomorphic three-form. In the superpotential W, the flux contribution to 



W (Ui, S) depends on the dilation S and the h ' complex structure moduli Ui (i — 1, 2, h 
while the non-perturbative terms for ft, 1 ' 1 Kahler moduli (A; = 1,2, .., h 1 ' 1 ) are introduced 
in [To]. The dependence of on Ui, S are suppressed. The model also includes the ol- 
correction (the £ term) to the Kahler potential 



11 



This model was originally proposed for 
the Large Volume Scenario [9] (see also 12 14]), and has been further analyzed in the search of 
de-Sitter vacua (oTlSl. 



Since we like to study the behavior of the expectation value of A when the number of 
complex structure moduli fields is large, we employ a simple model motivated by the orientifolded 
orbifolds of T 6 (7 , 15 , given by 



K, 



d+cs 



In (S + S) - ln {Ui + U) 



i=i 



h 2 ' 



1.2) 



W (Ui, S) =a + ^ biUi - S c 2 + diUi 



1=1 



1=1 



where 6j and di are (real) flux parameters that may be treated as random variables with 
smooth probability distributions that allow the zero values. Here we are interested in the physical 
A (instead of, say, the bare A), so the model should include all appropriate non-perturbative 
effects, a' corrections as well as radiative corrections. We see that the above simplied model 



(1.1) includes non-perturbative A k terms to stabilize the Kahler moduli and the a' correction 



£ term to lift the solution to de-Sitter space. In the same spirit, all parameters in the model, 



in particular the coupling parameters Cj, bi and di in Wo (1.2), should be treated as physical 
parameters that have included all relevant corrections. 



Following the analysis of 16 20 , we expect that the probability of an extremum at positive 



vacuum energy to be a classically stable solution will be Gaussianly suppressed (see also an 
estimation at supersymmetric AdS in SUGRA |21|). So the existence of solutions will put 
constraints on the parameters and we take s = Re(S) = l/g s > 1 for weak coupling. It 
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Figure 1: The probability distribution P(wq) for wo (1-3) with h 2 ' 1 — 1, 3, 5, 7, 9 number of complex structure 
moduli respectively, where wq is the value of the superpotential Wq at the supersymmetric solution. Here we 
present P(wq) in the range —10 < wq < (J_ 10 P{wo)dwo — 1) after solving for the complex structure moduli 
and the dilation at the supersymmetric point satisfying s > 1 for weak string coupling and Uj > for the Kiihler 
potential to be real. After some transient behaviors at small h 2,1 , the peaking (at wq — 0) of P(wq) strengthens 
as h 2,1 increases. In the last figure, we show the behavior of the expectation value (l^ol). 



turns out to be a good approximation to stabilize the Ui and S at a supersymmetric minimum 
(DsWq = Djj.Wq = 0) before turning on the corrections for the stabilization of the Kahler 
moduli, which then breaks supersymmetry. At the supersymmetric minimum with their axionic 
components sitting at zero, Ui 
are determined and 

w = W-olmin = - „^ 2 .i r. — r. r, (1-3) 



Re(Ui) > (required for fQ+cs (1.2) to stay real) and s > 1 
2(ci + sc 2 )Il£i 1 (l-sr i ) 



E?=i(i + ^)n^(i 



sr,- 



where r\ = di/bi and s is given as a function of the real random parameters c±, C2, rj. 

In general, the parameters will be fixed constants within the supergravity framework. How- 
ever, it is the flux compactification property in string theory that allows us to compare solutions 
with different choices of values for the parameters (e.g., different choices of fluxes may stabilize 
the complex structure moduli and dilaton differently). So we are justified to treat the parameters 
as variables with some suitable probability distributions and study the consequences. In this 
sense, the mechanism we suggest here can be considered as a stringy mechanism. As we sweep 
through the (e.g., uniform) distributions for the parameters c\,C2,bi,di, we see that each factor 
in the numerator of ( |1.3 ) easily passes through zero. So we expect the probability distribution 
P(wq) of the value for the superpotential w = Wn\mm to become more peaked at w Q — > as 
the number h 2 ' 1 of complex structure moduli increases (e.g., to hundreds). This behavior is 
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Figure 2: The probability distribution P(Wq) of Wq (1.2) when the real components of the dilation and the 



complex structure moduli are also treated as random variables. Here we assume uniform distributions with 
ranges -1 < c i: fy, di < 1, 1 < s < 5, and < Ui < 1. The first 2 figures are for P(W ) for ft, 2 ' 1 = 7, 9. The third 
figure is for P(D Ui W ) of D Vi W for ft 2 ' 1 = 7. 



illustrated in Figure [Tj 



It is interesting to compare P(wq) of wq (1.3) after stabilization with the probability distri- 



bution P(Wq) of Wq (1.2) before stabilization. Let us choose the same uniform distributions for 



q, 6j, G?j used in Figure [TJ but instead of solving for u,- t and s, we treat them as random variables 
with uniform distributions with range < Ui < 1 and 1 < s < 5 (while the axionic modes are 
suppressed). With these uniform distributions, we show P(Wo) in Figure [2j it is clearly smooth 
at Wo = 0. We note that the probability distributions P(DuiWo) for the variables Dy-Wo are 
also smooth with little or no preference for the zero value, as expected. If P(Wq) and P(Djj.Wo) 
are truly independent distributions, then fixing D v .Wo = (and D S W = 0) will not change 
the distribution P(W ); that is, P(w ) should be the same as P(Wq). However, this is clearly 
not the case when we compare Figure [l] (for P(wq)) and Figure [2] (for P(Wo)). This is because 
the two "random" distributions P{Wo) and P^Du^Wo) are actually correlated. In general, given 
any specific model, we expect correlations between these two distributions (as well as those for 
DDWq and DDDW ). It is this correlation that leads to the peaking feature in P(wo) in Figure 

m 

Next we insert the value wq into the superpotential W which is then inserted into the 
potential V. This allows us to find A for the de-Sitter meta-stable vacua via the stabilization 
of the Kahler moduli. This hierarchical setup effectively reduces the number of moduli fields 
when we reach the energy level for the Kahler moduli stabilization with determined s and the 
complex structure moduli itj, and thus helps to enhance the probability for dS vacua. Since the 
SUGRA approximation is valid only if the scale of the potential V is around or below the Planck 
scale, we have to restrict ourselves to such a valid set of potentials V in the determination of 
A. Among these meta-stable vacua, we expect that the probability distribution P(A) for A will 
also peak at A = 0. The numerical result on P(A) for the single Kahler modulus scenario is 
shown in Figure [3] 
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Figure 3: The probability distribution -P(A) of the cosmological constant A at meta-stable vacua as a function 
of h 2 ' 1 = 2,5,8 number of complex structure moduli and a single Kahler modulus (h ' = 1). Although the 
range is < A < 1, the probability distributions for only < A < 10 -3 are shown. P(A) becomes more peaked 
at A = as h 2 ' 1 increases. 



I. 
■ • ■ 



w 



• • t • a 



A 



(,2,1 



Figure 4: The LHS shows the average (\wq\) (red circles), |u>o| 80% (blue squares) and |wo| 10% (green diamonds) 
as a function of the number h 2 ' 1 of complex structure moduli. The RHS shows (A) (red circles), A 80% (blue 
squares) and A 10% (green diamonds) as a function of h 2 ' 1 . That is, there is a Y% probability that A of the 
meta-stable vacua will fall in the range A Y% > A > 0. For example, A 10% ~ 4.8 x 10~ 28 at h 2 ' 1 = 20 and 
~ 1.5 x 10 -41 at h 2 ' 1 — 30 (green diamonds). The values are for the case where the 6, parameters are fixed. 

Although the peaking property of -P(A) is essential for the preference of a very small A, it 
is not enough. We find that the expectation value (A) for de-Sitter vacua more or less stays 
constant (or changes little) as the number of complex structure moduli increases. That is, -P(A) 
in this scenario has a long tail outside the A ~ region. 

However, there are more than one way to implement the distributions for the random pa- 
rameters. Since Ui = Wo/2bi(l — sr^), allowing fej = implies that we are looking at vacua where 
some Ui — > oo, which may not be valid within our approximation. So, as an alternative, let us 
restrict the values of bi to a range that excludes zero, say > —f m i n > h > —f ma x ■ It turns 
out that the qualitative physics does not change if we simply restrict bi to a fixed value, say 
bi = — f < 0, where we may consider — / to be the mid-point of the range for bi. We shall choose 
/ so that 90% of the resulting A falls in the range 1 > A > 0, which we keep (our analysis is 
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Figure 5: The probability distribution P(A) for the ft 1,1 = 1, 3, 5 Kahler moduli cases with uniformly distributed 
Wo,Ai. Note that P(A = 0) is increasing (slowly) as h ' increases. 



probably not valid for larger A). In this scenario, both P(wo) and P(A) are peaked at zero, not 
too different from Figures [l|and§ while, for A > 0, (A) ~ e - 2 - 56/l2,1+7 - 40 for h 2 ' 1 < 10. However, 
the drop of (A) slows down appreciably as h 2,1 > 10 (see Figure [4]). These examples show that 
the behavior of (A) is very sensitive to the details of the input distributions of the parameters. 
On the other hand, the peaking (divergent or not) behavior of P(A) as A — > is quite robust. 

To get a better feeling of the peaking property of -P(A) for A > 0, let us introduce A Y% , 

t-A Y% 

defined by J P(A) dA = Y%. That is, there is a Y% probability that A of the meta-stable 

vacua will fall in the range A Y% > A > 0. (Note that A 100% ~ 1.) In Fi gure[4j we show (A), A 80% 

and A 10% as a function of the number of complex structure moduli. For example, at h 2 ' 1 = 20, 

while (A) ~ 1.2 x HT 8 , A 80% ~ 2.0 x 10^ 19 and A 10% ~ 4.8 x 10~ 28 . For h 2 ^ = 30, we find that 
A io% ^ 15 x 10 -4i. that igj there 

is a 10% chance that A < 1.5 x 10 41 in this case. If we take a 
cutoff lower than the 90% used, the resulting A values will be bigger (i.e., less impressive). For 
comparison, we also show in Figure [Z] the average (wq) and the corresponding w 80% and wl° % . 

We also consider the multi-Kahler moduli cases. In Figure |HJ we show P(A) for h 1 ' 1 = 1, 2, 3 
Kahler moduli cases (where the complex structure moduli and dilation are suppressed). We 
see that the peaking behavior also strengthens as the number of Kahler moduli increases. For 
positive A, we have (A) ~ 0.00110 (/i 1,1 ) -0 ' 226 e _0 - 0268/i1,1 , where uniformly distributed W , Ai 
are assumed. Here, the coefficient in the exponent is actually too small to be taken seriously. 
However, if we have sharply peaked distributions in Wo and A^ a clear exponential suppression 
in (A) as a function of h 1,1 is expected. We shall illustrate these features with some simplified 
models. 

In short, we see that P(A) becomes more peaked at A = as the number of moduli increases. 
This feature is robust. On the other hand, the behavior of (A) (or a/ (A 2 )) is sensitive to the 
way we implement the random parameter distributions and other inputs. The model studied 
here is chosen for its simplicity. The numerical data presented here are for the specific scenarios 
described. There may be other interesting scenarios within this simple model. To go beyond, 
it is important to examine more realistic models to check (1) whether the peaking of P(A) at 
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A = is as robust as we believe, and (2) what models will lead to an exponentially decreasing 
(|A|) (or A y% ) as a function of the number of moduli. A better understanding of the underlying 
physics should allow us to determine the value of (|A|) or A Y% as a function of the number of 
moduli. An exponentially decreasing (|A|) or A Y% may be used as a criteria to select particular 
regions of the landscape. 

In general, the probability that an extremum happens to be a meta-stable vacuum is Gaus- 
sianly unlikely as the number of moduli increases, as discussed in 16 21 . It is our goal to find 
scenarios where an exponentially small A for a meta-stable vacuum is preferred when the num- 
ber of complex structure moduli is large, as one expects to be the case in any realistic model. 
In fact, some of the commonly studied Calabi-Yau 3-folds with one Kahler modulus are smooth 
hypersurfaces in CP 4 with hundreds of complex structure moduli. So an exponentially small A 
hopefully may be naturally realized in this framework. 

Note that having a very small A does not necessarily solve the A problem completely. In 
this setup, the expectation value of the mass (squared) of the lightest modulus tends to decrease 
more quickly than (A). So we may have the cosmological moduli problem 22 24 . It remains 
to be seen whether the light modulus mass is really a phenomenological problem as one may 
intuitively envision. 

In summary, we find that the study of multi-moduli scenarios is a promising direction in 
the search for the reason why the observed A is so small. Tremendous amounts of effort have 
been spent on the search for a string vacuum with a standard model of three families of quarks 
and leptons. We believe it is worthwhile to search for regions of vacua in the stringy cosmic 
landscape that have a chance of providing a naturally small positive A. Combining these two 
criteria in a search may be more effective than the searches carried out so far. On a separate 
note, we hope that this product distribution idea can find applications in other phenomena. 

The rest of the paper is organized as follows. Section [2] reviews some simple probabilistic 
properties of functions of random variables and the single Kahler modulus model. Section [3] 
discusses the multi-complex structure moduli case. Here we find the supersymmetric solutions 
for the complex structure moduli and the dilation and then insert them into the single Kahler 
modulus model to solve for A. Section [4] discusses the multi- Kahler cases. Here, we also consider 
some simplified versions of the general multi- moduli cases. Section [5] presents the summary and 
discussions. Some details are relegated to the appendices. 



2 Background and Review 

We begin with a review of the basic probability distribution of a function of random variables 
when the distributions of the random variables are given. In particular, we emphasize on the 
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peaking property in the resulting probability distribution, since this is relatively insensitive to 
the details. We explain how a product of random variables gets a sharply peaked (at zero) 
probability distribution. We also point out the importance of the overall distributions to the 
expectation values. Then we review how this feature emerges for the probability distribution of 
A in the single Kahler modulus model in a class of models in Type IIB string theory well studied 
already, basically following |3|. We shall use this model to provide the framework to study the 
multiple complex structure moduli scenarios. In appendix |A| we give some field theory toy 
models to illustrate some of the probabilistic properties that one may encounter in this analysis. 



2.1 Probability distribution of functions of random variables 

Suppose we have n random variables yi (i = 1,2, ••• ,n), each with probability distribution 
Pi(yi), where f Pi{xji)dyi = 1. Let us consider an arbitrary function z = f(yi,yz, 4 4 4 ,Vn)- Then 
the probability distribution P(z) of z is given by 

P(z) = / dy 1 dy 2 ■ ■ ■ dy n P 1 (y 1 )P 2 (y 2 ) ■ ■ ■ P n {y n ) &{f{Vi) ~ A- 

(2.1) 



J P{z)dz =1 



so the probability distribution P(z) of z can always be properly normalized, even when P(z) 
diverges at z = and/or elsewhere. 

Consider the following simple example: z = y\y 2 • • -y n an d each random parameter y^ obeys 
the uniform distribution P(yi) — 1/L with range < yi < L; then the probability distribution 
P{z) of z is given by (for < z < L n ) 

Therefore P(z) of a product of random parameters has a mild divergent peak at z — 0. This 
divergent behavior strengthens as n increases. Actually this logarithmic divergent behavior of 



P(z) at z = (2.2) is present for any yi distribution that smoothly includes yi = 0. 

Given the probability distribution for yi, we can determine its expectation value. Here, for 
z = yiV2 4 4 4 yn, the expectation value of the pth moment (z p ) of z is given by 

(z p ) = {vl) ■ • ■ (y p n ) (2.3) 

following from the property of Mellin integral transformation. Since (yi) — L/2 for the uniform 
distribution, we have 

( z )= (-X = e -"M2/£). (2.4) 
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For L < 2, we see that (z) decreases exponentially as n increases; so the divergent behavior at 
z = is correlated to the smallness of (z). However, (z) increases exponentially if L > 2. In 
this case, the long tail (at large z) of P(z) drives (z). In short, we see that the behavior of (z) 
(as a function of n) is very sensitive to the property of P{z), especially at large z. 

To conclude, we like to mention that P{z) is in general smooth at z = if z is a sum of terms 
that involve independent random variables (see e.g. |3| for more details). The problem now is 
to find the functional form of the cosmological constant A of meta-stable vacua in terms of the 
random parameters entering the model and see whether its distribution has the nice peaking 
property at zero and whether its expectation value or A Y% decreases exponentially fast as a 
function of the number of moduli. 



2.2 The single Kahler modulus model 

We briefly review the analytic study of the probability distribution P(A) in the single Kahler 
modulus scenario in type IIB string theory |3j|. Although no potential is generated for the 
Kahler modulus at the classical tree level, non-perturbative terms can be introduced into the 
superpotential W and a'-corrections into the Kahler potential to break the no-scale structure, so 
that a non-trivial potential is generated at next leading order. The tree level potential basically 
controls the stabilization of complex structure moduli and dilaton. Here we simply assume the 
complex structure sector is stabilized at high-energy levels and therefore is believed to have 
been integrated out in the low energy effective model. Later, we shall introduce the complex 
structure sector, solve it at the supersymmetric point and insert the solution into this model to 
solve for the stabilization of the Kahler modulus which then breaks supersymmetry. 

The model of interest is given in [5148] , also known as Kahler uplifting model : 



(2.5) 



K = - 2 In 

W =W + A ie ~ aiTl 



( V + f)' V = 7 i(Ti + f 1 ) 3 / 2 , 



where the a'-correction related to the £ term [11] and non-perturbative correction with the 
parameter A\ llO] are given. The parameters Wo, A\ will be treated as random variables. We 
work within the approximation 



— 1 
V ' 



-a x T x 



W 



« 1, (2-6) 



to simplify the analytical study [7]. We shall see that this approximation is quite reasonable 
for the analysis of small A. The first approximation is a sort of the large volume approximation 
which is powerful for suppressing the other unwanted stringy corrections. 
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Figure 6: The potential V for different values of C. Only the bracket part in V (2.7) is shown. The V H defined 
at C ~ 3.89 gives roughly the height of the potential barrier. 



Since the imaginary part of the Kahler modulus T\ has a cosine type of potential, the extremal 
condition for this direction is satisfied when ImTx = 0. Therefore we focus only on real part 
fi = ReTi. 

Now the potential simplifies to [FI 



V 



c 



■W Q a\A x ( 2C e 



-27W Q £af 2 
64v^7i^i 



9zf 



x 



351 



(2.7) 



Xi = aiti 



The stability condition d% V > at the extrema = with respect to x\ is easy to analyze, 
and we get the parameter range for stable positive A: 



C Q <C<C X -+ 3.65 < C < 3.89, 



(2.{ 



where the lower bound is given by positivity of the minimum of V, while the upper bound 
is given by the stability constraint. If we satisfy the combination of parameters C inside this 
region, there is a stable solution in the range 2.50 < x\ < 3.11 at A > 0. Up to an overall factor, 



the potential V (2.7) is shown in Figures 



Our interest is to study the statistical property at small A. In this approximation, the 
extremal condition gives us the linear relation between x\ and C by 



Xi 



-5/2 



e>' 2 C - 2. 



(2.9) 
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Figure 7: The comparison between the analytic result ( 2.12 ) (solid line) and the numerical result of the potential 



(2.51 (histogram). Here we set a± = 71 = £ = 1 and assume uniform distribution for parameters W$,A-\_. The 
plot in LHS shows the probability distribution P(A) of A, while the plot in RHS is the probability distribution 
P(ln A) vs In A. They are related via P(A) dA = P(ln A) rfln A. 



and we have 



A=n 



9/2 



-W Q a\Ai 



(C - 3.65) . 



(2.10) 



Looking at V (2.10), we note that the height of the potential barrier is relatively unchanged 



as A varies from negative values to positive values, as shown in Figure [6] So we may use the 
barrier height as a measure of the scale of the potential V, 

9/2 



Vj 



H 



(0.24). 



(2.11) 



In general, the SUGRA approximation is valid if the scale of V is less than the Planck (or string) 
scale. We shall require Vh < 1 as the condition for the validity of our approximation. 

In finding the probability distribution P(A), the peaking behavior of the distribution comes 
essentially from this Vh- Since we would like to focus on the shape of distribution, we set the 
maximal value to be Vff| max = 1 here, such that the effective theory is safely applicable under 
the Planck scale. For further simplification, we may set a 1; 71, £ to constants and see the effect of 
two random variables W , Ai on P(A). (In |3|, we also discuss cases where the other parameters 
are treated as random variables.) 

If we simply assume < — Wo, A\ < 1 obeying uniform distribution, the probability distri- 
bution can be obtained (the detail is relegated to appendix [B| , 



P(A) 



2500^5 
3SC 



In 



5C 



2500\/5 A 



(2.12) 



where 5C = {C\ — Cq)/C\ = 0.0617. So we clearly see the divergent behavior as A — > 0. 
The stabilization dynamics is responsible for the divergent behavior in -P(A). As we state in 
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the introduction, the peaking behavior of P(A) is robust. Note that this divergent behavior 
is present for negative A as well. For A < 0, there is another branch of supersymmetric AdS 
solutions; so .P(A) for A < is not a reflection of the P(A) shown above. 



Let us compare P(A) obtained analytically in the approximate potential (2.10) with that 
obtained numerically for the potential (2.5), as shown in Figure [7J From the two figures there, 
we see their nice agreement with each other, especially at small A. 



3 Multiple complex structure moduli cases 

So far we have reviewed the case with a single Kahler modulus assuming moduli stabilization 
of the complex structure and the dilaton at higher energy scale. In the actual models, we 
have the dilaton S and many complex structure moduli Ui (i = 1, 2, h 2 ' 1 = N cs ) which are 
assumed to be stabilized at a higher scale with fluxes as sources. Here, we analyze in more 
detail their stabilization and see how it affects the probability distribution of Wq. In (25], a 
random distribution for the superpotential (and its derivatives) is assumed to investigate the 
entire landscape of SUGRA. Here, we probe parts of the SUGRA landscape in greater details. 

Throughout this section, we focus on models of the form: 

K cs+d = - In (S + S) - ^ In (Ui + ft) , 

(3.1) 

S. 



h 2 ' 1 / h 2 ' 1 \ 

W =Ci +J2 b iUi~ [C2 + J2 d jUj 
i=l \ i=l / 



where Cj, b{, di are parameters from fluxes. This model is motivated by the orientifolded orbifolds 
of T 6 (7, 15 , and is used for the purpose to see the tendency of the probability distribution of 



Wq as a function of h > here. Supersymmetry breaking comes from the interplay between the 



a' correction £ term and the non-perturbative term reviewed in section 2^ So it is a good 
approximation to solve for Ui and S at the supersymmetric point and then insert their values 
into the full system to solve for the stabilization of the Kahler moduli. It is worth noting 



that we basically require £ > (1.1) for the stability analyses of Kahler moduli later, meaning 
that the number of complex structure modulii is greater than the number of Kahler moduli: 
x(M) = 2(h 1 ' 1 — h 2 ' 1 ) < 0. As pointed out in [7], this approximation is valid when the volume is 
large, as we always assume, and when h 2 ' 1 ^> 1, which is precisely the region we are interested 
in. 
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3.1 Super symmetric Minimum 



Following the approach in (7], we stabilize the dilation S = s + ia and complex structure moduli 
Uj = Uj at the SUSY point to find Wolmin and then insert it into the superpotential W for 
the multi-Kahler moduli stabilization where SUSY will be broken due to the non-perturbative 
effects. To find the supersymmetric values for Ui and S, we have DgWo = dsWo + dsKWo = 
and Du.Wq = 0, or 



Re(ZWo) = Re(D Ui W ) = 0, 
Im(ZWo) =lm(D Ui W ) = 0. 



(3.2) 
(3.3) 



For s > 1, (3.2), (3.3) have solutions so the imaginary modes decouple in (3.2). So it is most 
convenient to work with the a = = solution for the imaginary modes and focus on the real 
parts of the moduli fields in the region s > 1 (for weak coupling) and Ui > (to ensure that 
K cs+< x stays real). Now Wq becomes (defining N cs = h 2,1 for convenience) 



Ac 



Wn 



(ci - sc 2 ) + /Xh - sdi)Ui 



(3.4) 



and D S W = and Djj.Wq = yield (assuming s ^ and ^ 0) 



iVr, 



(ci + sc 2 ) + y~)(bj + sdi)ui = 0, 



JVc; 



(ci - sc 2 ) - (hi - sdi)m + /;(bj - sdj)uj = 0. 



(3.5) 
(3.6) 



Now A^cs = and N cs = 1 are special cases and are easy to solve. So let N cs > 2. We shall 
see that iV cs = 2 is also special. (3.6) immediately gives 



v = (&i — sdx)ui = (b 2 - sd 2 )u 2 



(b Ncs - sd Ncs )u Ncs 



(3.7) 



So Ui are solved in terms of s and one of them, say or equivalently, v. Then (3.6) becomes 

(ci - sc 2 ) + (iV cs - 2)v = 0, (3.8) 

and thus 

Ci — SC2 

W^oLmin = -2^T^| = (3 ' 9) 

where the second form is also applicable to the iV cs ^ cases. We see that s = c\/c 2 for iV cs = 2. 



Now we have two equations (3.5) and (3.8) for two unknowns, v (or u\) and s, and we have 
for N cs ^ 2, 



(N cs - 2) 



Ci + sc 2 
Ci - sc 2 



Ebi + sdi 1 + sr. 

h- - S d; ~ 



i=l 



i=l 



ST, 



(3.10) 
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where r% = di/bi. In general, we cannot solve analytically for s. However, once s is determined 



by (3.10), we can express Wolmm^cs) m terms of s, 

«» s w.|-,(Ag = - N 2(C1 + SC2) = _ 2(c 1 + sc 2 )ut 1 (b,-sd,) 

E£; a + <"\)/(l - sr.) ESft + «4)nj^(6j - «*,) 

This is the formula quoted in the Introduction. Recall that we like s to be large and positive 
say s « O(10). For the iV cs = 0, 1, 2 cases: 

^ |min(iV cs = 0) =2ci, s = --, 

c 2 



ty.U(JV.= l)=-2 ^±^^ , S = (3.12) 

Oi + SCti V C20ti 

w i c /v _ ^ _ (ci + sc2)(6i-sdi) (62-^2) o _ Ci 

1^0|min(iV cs - 2) — — — , S - — . 

°1°2 — s «1U2 c 2 

In appendix [Cj we further discuss the solution for s in special cases. 

Treating Cj, bj and dj as random parameters, we allow them to take values such that each 
solution satisfies s > 1 and Ui > 0. We see that the each factor (bj — sdj) in the numerator 



in (3.11) is allowed to pass through zero. We will check numerically that there is actually a 
peaking behavior as the number of complex structure moduli increases. 

Before introducing the next-leading order terms (both non-perturbative and the £ term) to 
generate a potential for the Kahler moduli, the model has a no-scale structure so I-Dt^^I 2 = 
3 1 Vt^ 1 2 . Since V = e Y2i \DiW\ 2 (with i running over complex structure moduli and dilaton) is 
semi-positive definite (with the supersymmetric minimum at V = 0), the Hessian (mass squared 
matrix) is semi-positively definite, owing to the no-scale structure. On the other hand, some 
of masses for the axionic modes turn out to be zero as discussed in [7]. Presumably, the non- 
perturbative corrections will generate cosine-like potentials, thus stabilizing the axionic fields. 

The above values for s and u\ are solved at the supersymmetric point. Their values will be 
shifted after supersymmetry breaking [7]. The corrections are estimated of order 



5(s,Ui 



S,Ui 




(3.13) 



As we reviewed in section 



2.2 



the expansion of 0(£/V) is compatible with the analysis for smaller 
A, and thus we can neglect the small correction to the complex structure sector in our analysis. 
It is worth commenting about the situation that total summed over corrections with large N cs is 
not negligible. In this case, the extremal condition in the perturbation d^d^V ■ 5(f) j = —d^.SV 
suggests roughly 

6( s , Ui )<xO (-?-). (3.14) 



Therefore the correction may be suppressed by large N cs instead 
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3.2 Some properties of wq = Wo|min 



We see that w (3.11) together with the solution for s (3.10) have some interesting properties. 



Since c±,C2, bj, dj are random variables, so, for any fixed j, (bj — sdj) — » is surely allowed. 



(3.10) for s indicates that when (bj — sdj) — > for any j, the RHS blows up (in the absence 



of delicate cancellations) so LHS requires s — > C\jc2- Now, in the vicinity of bj 
s becomes independent of the other values of (b iy di) for i ^ j. 



sdj = 0, 



The factors (bi — sdi) appear in the numerator of w (3.11). Knowing that s — > Ci/c 2 
is insensitive to (bi, di) for i ^ j in the vicinity of (bj — sdj) = 0, we can work out the 
probability distribution P(z) for z = bi — sdi for i ^ j, 

P(z) = I I dhddi S ((bi - sdi) - z) . (3.15) 



o 



A careful treatment of the limits of integration yields, 



P(z) 




< z < 1, 

1 - s < z < 0, 
— s < z < 1 — s. 



(3.16) 



and P(z) = outside. This distribution has the shape of a trapezoid and P(0) = 1/s. (It 
follows that (z) = (1 — s)/2.) This implies that when (bj — sdj) — > 0, the probability that 
(bi — sdi) — > for any i ^ j is not suppressed. Furthermore, the individual (bi — sdi) — ^ 
approach are more or less independent, so the n = iV cs th order zero is within the parameter 
ranges we are interested in and we may intuitively expect a peaking behavior 



P(lU ) OC (hi | Wo |) ^ 1 as Wq — > 0. 

Let Ti = di/bi. At the symmetry point where all rj are equal, T\ 
see that W^olmm^cs) ( |3. 1 1 ) collapses to 

2(ci + sc 2 )(l 



r-2 



rN cl 



W | min (iV cl 



sr) 



(3.17) 



r, we 



(3.18) 



N cs (l + sr) 

As r sweeps pass 1/s, we obtain only a single zero in wq. For this symmetry point, we 



expect P(wq) to be smooth at wq = 0. That is, the peaking behavior (3.17) will disappear 



if we restrict ourselves to the symmetry point. So if the peaking behavior in P(wq) is to 
survive, it must come from the vicinity of the symmetry point, which has a bigger region of 
parameter space than that of the symmetry point. Away from this symmetry point, some 



zeros in the numerator of wq (3.11 ) may be canceled by coincident zeros in the denominator 



of Wq (3.11), depending on the choice of the random parameters. This may lead to some 



transient behavior in P(wq). Solution for s at the symmetry point is discussed in appendix 
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Constraint among the parameters comes into play as we limit the range of s so the weak 
coupling approximation is valid. This will also introduce some transient behavior in P(wq). 
This transient behavior should fade away as we go to large N cs . 



Note that both the expression for s (3.10) and the expression for u> (3.11) are invariant 



under the simultaneous rescaling 6j — > AA and di — > Xidi, so P(wq) should be unchanged 



when we change the range of both hi and di by a factor of Aj. Recall the example (2.3) 
where the expectation value strongly depends on the range of the random variables. This 
would be very unsatisfactory here since we do not know the reasonable ranges for (bi,di). 
Fortunately, because of this Aj scaling property, we see that both P(w ) and (wq) are 
insensitive to the rescaling of the (bi, di) ranges. We shall find some reasonable choices of 
ranges for hi when we come to the potential V and -P(A). Here, we restrict ourselves to 
Wq < so we simply need to look at only (w ) (instead of (wq)). Suppose we take the 
range of q to be [-r cutoff , +r cutoff ], then (w ) /r cutofr is independent of r cutoff and the ranges 
of bi, di. For example, for N cs = 5, we have (wo) /r cuto ff — —2. 

Overall, we expect that P(wq) peaks more sharply at wq = as N cs increases; however, we 
do not know analytically the precise peaking behavior of P(wq). 

So we expect that the expectation value (u> ) decreases as N cs increases. Because of the 



nontrivial forms of s (3.10) and Wq (3.11), it is crucial to find their behaviors numerically. 



3.3 Numerical comparison of wq 



In the above analysis, we see the product type of structure in Wolmin (3.11). However, since the 
parameters in the denominator are not independent of those in the numerator, it is important 
to perform numerical analysis and check whether the distribution does become more sharply 
peaked as we increase the number of complex structure moduli. 



Now we impose a randomness in the model (3.4). Given random values for ci, C2, di obeying 



uniform distribution, we solve supersymmetric conditions D$Wq = 0, Du.Wq = for real part of 
moduli fields, s,Ui. Here we restrict ourselves in —1 < c\,C2,bi,di < 1 for just a choice. We are 
also interested in the region where string coupling g s = = s~ l can be treated perturbatively, 
therefore we further impose s > 1 together with Ui > (which is required for a real-valued 
Kahler potential in our definition). 

Recall that the WqAi < condition is required for the Kahler moduli stabilization with 
positive A as we see in section 2.2| Because of the reflection symmetry, we may focus on the 



wq < region without losing generality here. It is worth commenting that all of the numerical 
solutions obtained just by solving the supersymmetric condition actually satisfy the positivity 
of the Hessian. Thus we can insert the solution in the model for the Kahler moduli stabilization. 
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Figure 8: The probability distribution P(w Q ) for -10 < w < 2 with fixed bi = -1 at iV cs = 1,3,5,7 
respectively. In the last figure, we show the behavior of the expectation value (|u>ol) f° r ^cs < 10. 



Inserting the solutions for s (3.10) and Ui ( |3.7[ ) into wo = Wo|mm (3.11), we obtain the 
distribution P(u>o) and its behavior. In Figure [TJ we show the resultant probability distribution 
P(wq) for various h 2,1 = N cs , where the dilaton dependence has been taken into account. Because 



of the non-trivial dependence in the denominator of (3.11), the peak at Wq = is smoothed out 



around iV cs = 3,4, but emerges at N cs > 5. For iV cs > 5, the peak at wq = becomes sharper 



as we increase N cs . Although Wq (3.4) takes the form of a sum of terms and therefore one 
might expect a sum distribution (which is smooth at Wq = 0, as shown in Figure [2]), instead 
of a product distribution (which typicaly peaks at 0), the data clearly shows a sharper peak 
in P(wo) (Figure [T| as we increase the number of complex structure moduli. This is nothing 
but the outcome of the correlation of the terms through the stabilization of complex structure 
moduli and dilaton. 

Although the distribution feature for the parameters ci, C2, bi, dt are reasonably well-motivated 
from the flux compactification because these parameters are from the 3- form fluxes H 3 ,F 3 , we 
like to consider other possibilities. Since Ui = w /2bi(l — sTi), allowing 6^ = implies that 
we are looking at vacua where some Ui — > oo, which may not be valid within our approxima- 
tion. So, as an alternative, let us restrict the values of 6« to a range that excludes zero, say 
> —fmin > b{ > —fm ax . It turns out that the qualitative physics does not change if we simply 
restrict bi to a fixed value, say bi = — f < 0, where we may consider — / to be the average value 
of bi with a range. Note that Ui stays finite as (1 — sri) — > 0, since wq — > in this limit as well. 

Again, we can solve for s > 1 and Ui > and determine wq. In Figure |8j we show the 
distributions of P(wq) with fixed bi = —1 while allowing the other parameters uniformly dis- 
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tributed with the range —1 < Cx,C2,d{ < 1. We again impose the weak string coupling constraint 
s > 1 as well as Mj > 0. As apparent from the comparison between Figure [T] and Figure [8} the 
distribution of w in this case is more peaked toward w = 0. Note that the distribution for 
wq > is quite different from that for wo < in this setup. This is because the sign fixing of b{ 
violates the reflection symmetry of wo (realized by flipping signs of all the parameters q, bi, di 
simultaneously in the model). 

We also estimate the expectation value (\w \) as a function of N cs , especially for N cs = 4 — 10. 
In an attempt to extract the large N cs behavior, we neglect the data points for iV cs < 4, because 
of the transient behavior at lower N cs . The data points (for N cs > 4) in Figure [8] are well 
fitted by (\wo\) = 60.6 e ~ om3Ncs . However, numerical data indicates that the drop of (|tuo|) as a 
function of N cs slows down as we go to N cs > 10 (see LHS of Figure [5]). 

The peaking behavior of P(|tOo|) rnay be quantified via a comparison of the values at a 
fixed percentage of \wq\ accumulating from the zero point. Let us introduce |u>o| y% defined 

I \Y% 

by lo P(\w \)d\w \ = Y%. That is, there is a Y% chance that \wq\ will be smaller than 
|iuo| Y% . In LHS of Figure [ij we show (|iwo|), |wo| 80% and |wo| 10% as a function of A^ cs in the range 
1 < N cs < 30. We see that, for iV cs > 3, they all decrease as the number of complex structure 
moduli increases. 



3.4 Inserting the wq solution into the single Kahler modulus model 

Now we consider the moduli stabilization in the presence of single Kahler modulus, given the 
data of wo, s, Ui obtained by solving supersymmetric conditions for complex structure and dilaton 
moduli. The potential we have in mind here is given by 

K = - 2 In + |j - ln(2s) - ^ ln(2 Ui ), 

2 C(3)\ 3/2o ,aq ( 3 - 19 ) 
3(2vr) E 

W =w + A ie ~ aiTl . 

We consider the stabilization of single Kahler moduli with the potential above, under the as- 
sumption of complex structure moduli stabilization at higher energy, which is not affected by the 
next leading order corrections for Kahler modulus stabilization. Again, we choose the solution 
where ImTi = and focus on t\ = ReTi. The potential above has additional contribution in 



V =7i(7i + Ti) 3/2 , i = - ,3/2 2^1,! _ h 2,^ 



the overall factor in addition to (2.5), which is estimated by 

-WqclIAx ( 2C e~ Xl 



V 



64^71^1 



,L i 



(3.20) 



C — , X\ — (X\t\, 
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where we approximate the potential up to linear order of £/V,v4ie Xi /Wq, and iV cs = h 2 ' 1 just 
for convenience. Here, we set a\ — 71 = 1 for simplicity. 

In this section, we present the data obtained for the following 4 cases : 

1. -1 < a, c 2 < 1, -fi < k, di < f h 

2. -1 < ci, c 2 , Ti = bi/di < 1 - f 2 < bi < f 2 , 

3. -1 < ci, c 2 < 1, k = ±/ 3 , -/ 3 < dj < / 3 , 

4. -1 < ci, c 2 < 1, k = -f 4 , -fi <di< / 4 , 

with uniform distributions for the ranges of the random parameters as shown. Note that we 
impose uniform distribution for r« instead of di in Case 2. The positively defined parameter fj 



will be fixed later. For each case, we insert the given distributions into (3.1) and solve for the 



supersymmetric conditions (3.2) for s,Ui,w , and collect sets of data for s > 1, Ui > 0. The 



invariance under rescaling b^di simultaneously in (3.11), the peaking behavior of w < as well 
as the distribution itself is independent of fj. 

As we have discussed in the previous subsection, we observe the peaking behavior in wq as 
iV cs increases. However, we also see that P(ui) are also peaked toward Ui = 0, and therefore the 



coefficient of the potential V{3.20) goes essentially infinity, meaning that the resultant A itself 
may diverge. Since we use the effective potential for the 4D gravity analysis, the maximal value 
of the potential should be within the Planck scale for it to be valid within the approximation 
used. Since V(xi = 0) blows up, it makes more sense to use the barrier height as a measure of 
the scale of the potential. The easiest way to implement this constraint is for the barrier height 



(2.11) to be no higher than the Planck scale, or 



V H = I ( 2 -X /2 (0.24) <1. (3.21) 

That is, we keep only those Vh that are less than unity. Here we give A\ an uniform distribution 
in the range — 1 < A\ < 1. To enhance the data sample without changing the underlying physics, 
we may simply restrict ourselves to small values of fi, where Vh oc /j^ 8 - There may be a number 
of different ways to choose the cutoff point. Here we employ a specific way to achieve this: we 
count the number of data starting from Vh = 0, and we choose fi so that 90% of Vh < 1, then 
we set this as the cutoff point and keep only this Vh < 1- This implies that we restrict ourselves 
to the cases where A < 1. Note that we do not lose generality by the restriction wq < for 
Case 1-3, since we have the reflection symmetry of wq by changing the signs of all parameters 
Ci,bi,di, which remains solutions for s, Ui unchanged. Changing the 90% cutoff to a higher or 
slightly lower value does not change the qualitative picture. 

Let us estimate fj in each case. To keep 90% of Vh < 1, we need to choose a smaller fj as 
A^ cs increases. In LHS of Figure [9j we show the plot for the value of fj required to normalize the 
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Figure 9: The estimated values of fk and (A) in Case 1 (purple, rectangular) and Case 4 (blue, circle). The 
fitting curves here are f x ~ e -5.33+i6.9/w ca ( red Qn LHS ^ j 4 _ e -s.i6+25.7/JV cs ( blue Qn LH S) and (A) 4 ~ 
e -2.56Ar cs +7.40 (bi ue on RHS). However, the drop of (A) as a function of AT CS slows down appreciably for larger 
values of Ar cs . 

cutoff point of Vh in Case 1, 4 just for illustration. For N cs > 4, fj may be estimated to be : 

f l „ e -5.33+16.9/iV cs ^ ^ ^ e -4.62+20.5/JV cs ^ ^ ^ e -5.78+19.6/JV„. ^ ^ e -8.16+25.7/iV CB _ fa £2) 

We see that /j C 1 at larger N cs . 



Now we are ready to plug the data into (3.20). Since we consider just the single Kahler 



modulus stabilization with the sets of data inputs, the stabilization mechanism follows that 



reviewed in section 2.2, where the combined parameters obey 3.65 < C < 3.89 for metastable 
dS vacua. First we give A\ an uniform distribution with range —1 < A\ < 1. Next we choose the 
data for one of the four cases mentioned above for each fixed N cs . The probability distribution 
-P(A) of A in Case 1 is illustrated in Figure [3j We see that there is the clear peaking behavior 
toward A = 0, which is getting sharper as A^ cs increases. However, we do not see any clear trend 
for (A) in RHS of Figure [9] This is probably because there is a non-trivial long tail in large A 
in -P(A) which contribute significantly to (A). In Case 2 and 3, we see a clear peaking of P(A) 
at A = 0, similar to that in Case 1; but again (A) behaves quite similarly to that in Case 1, 
without a clear trend one way or another. 

Next, let us consider Case 4, with b{ = —f± and wo «s 0. The result is shown in Figure [9j The 
peaking behavior of P(A) at A = is again evident, while (A) has a rather different behavior 
than the other 3 cases. The (A) in Case 4 is shown in Figure |4j We find that it is roughly given 
by 

(A) 4 ~ e - 2 - 567Vcs+7 - 40 (3.23) 

for A^ cs < 10. However, the drop in (A) slows appreciably for A^ cs > 10 as showed in RHS of 
Figure|4j Still, this example illustrates the sensitivity of (A) to the input ranges and distributions 
of the random parameters. 

Similarly to Wq, we can quantify the peaking behavior of P(A) by estimating the value at a 
fixed percentage of data counting from A = 0. Let us introduce A , defined by J„ -P(A) dA = 
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Y%. That is, there is a Y% chance that A will fall in the range A Y% > A > 0. In Figure [4J we 
show (A), A 80% and A 10% as a function of the number of complex structure moduli. For example, 
at N cs (= h 2 ' 1 ) = 20, while (A) ~ 1.2 x 10" 8 , we have A 80% ~ 2.0 x 10~ 19 and A 10% ~ 4.8 x 10~ 28 . 
Although the drop of (A) is quite mild after N cs > 10, we see that A 80% and A 10% continue to 
drop appreciably. This means that most A have small values, while there exist a long tail A < 1 
in the distribution P(A). 

By choosing smaller allows us to accumulate more data. What happens if we impose a 
different cutoff on Vjj- As an example, let us choose 6j so that, instead of 90%, only 20% of the 
potentials have Vh < 1. That is, the potentials are less squeezed to small values in the 20% 
cutoff scenario. In this case, we find that, for A" cs = 20, A 80% ~ 8.6 x 10~ 6 and A 10% ~ 1.2 x 10" 9 , 
which are substantially bigger than those for the 90% cutoff case. Clearly we have to understand 
this squeezing issue better. 

So far we have considered the physical quantities A and the related potential barrier Vh in 
this subsection. However, one may also consider, instead of Vh, the quantity e^^l Wo| 2 > which 
is invariant under the Kahler transformation. The Kahler invariant quantity at a minimum is 
given by e is:cs + d |ll / o| 2 = WQ/(2 Nca+1 sY[v>i)- Since this quantity may diverge due to the Ui in the 
denominator (as in the case of Vh), we need to introduce the cutoff to make this quantity within 
Planck scale (this quantity is a term in the potential). If we follow the same procedure as that 
for Vh, the distribution of the Kahler invariant quantity at each N cs suggests sharply peaked 
distribution, which is qualitatively similar to what we have seen already for P(A). 



3.5 Some remarks 



A few remarks may be in order here: 



Note that the "product" form emerges in the superpotential value w (3.11), despite the 



fact that the parameters 6, and di appear in the "sum" form in the original superpotential 



(3.1 ). As pointed out in 3], any linear combination of a set of random variables, as naively 



the case in (3.1), will not yield a peaking behavior. That is, Wo will be smooth at zero if 



the Ui, the real part of Ui in Wq are treated as independent arbitrary constants. When we 
go to the supersymmetric point for the minimum, m become functions of 6j and di, which 



leads to the product form in (3.11). It is the (complex structure and dilation) moduli 
stabilization dynamics that converts the "sum" form ( 3.1[ ) to the "product" form (3.11), 
which is the key to the peaking property of P(wq) and -P(A), the latter of which is necessary 
for the stringy mechanism for a naturally small A to work. 

The actual values of Ui and s will surely be shifted away from the supersymmetric point 
after supersymmetry breaking. The correction of back reaction to the set of solutions 
is suppressed due to the suppression of large volume and larger Wq compared to the non- 
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perturbative terms if we work in this parameter region [7] . Even if we include the corrections 
in Wo, the modification would be negligible and therefore the feature of peaked behaviors 
should remain intact. 



It is possible (even likely) that the coefficient in (1.1) emerges as a consequence of the 
stabilization of higher scale moduli. In this case, we expect the probability distribution 
P(Ak) to be peaked at = 0. This will enhance the peakiness of P (A), so fewer complex 
structure moduli may be needed for a small enough A. 

Both the model and its solution we study here are simple enough for a clear analysis. 
This opens some questions on the simplifications we are taking. There are additional 
interaction terms in both W and K coming from higher order corrections. Based on the 
limited experience we have developed so far, one may hope that those (higher a' or stringy 
loop 13 26,27]) corrections will lead to tighter interactions among the moduli and so will 
actually strengthen the peaking behaviors. However, this does not necessarily imply an 
exponentially small A. 



4 Multi-Kahler moduli cases 

A model for multiple Kahler moduli stabilization at positive A is even suggested especially with 
non-perturbative effects in superpotential and a' correction in Kahler, considered in [5}j8] or 
known as Kahler uplifting model. The model basically assumes the stabilization of the complex 
structure moduli and dilaton-axion moduli at high scale, while the Kahler moduli stabilization 
is achieved at low energy scale with the non-perturbative effect and the a'-correction. Together 
with the non-trivial constant Wq in the superpotential, given after complex and dilaton moduli 
stabilization, we can have a small positive minimum as a result of Kahler moduli stabilization. 

We consider the supergravity effective potential which is obtained after the Calabi-Yau com- 
pactification in type IIB: 

N K 

K = - 2 In ( V + |) , V = 7i(7i + T^ 2 - ^ 7^ + T^ 2 , 

i=2 



w 




(4.1) 



where we have again introduced the non-perturbative terms proportional to Ai for all of Nk = 
h 1 ' 1 Kahler moduli, and the a' correction defined £ = — ((3)x(M)/ (4\/2(2it) 3 ). We basically 
require £ > for the stability analyses later, meaning that the number of complex structure 
moduli is greater than the number of Kahler moduli: x{M) = 2(h 1,1 — h 2,1 ) < 0. Wo comes from 
the flux contributions, and generically does not depend on Kahler moduli (no-scale structure). 
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We assume that the non-perturbative effect and the a' correction are small enough compared 



to the other quantities as well as in section 2.2 so that we analyze the potential up to linear 
orders of 



Ate' 



< 1. (4.2) 



W 

These assumptions are quite natural since the a' correction can be treated to be small in per- 
turbation theory. As we see in the previous section, this perturbative analysis agrees with the 
analysis for small A. There are essentially two important outcomes of these assumptions: (1) 
suppression of the off-diagonal components between the Kahler components and the complex 
structure/dilaton components in the mass matrix, and (2) simplification of the potential for the 
Kahler moduli since we can neglect the next-to-next leading order in the potential with respect 



to the approximation (4.2) 



We focus on the real part of Kahler moduli Re T = tj because the imaginary part always 
has a potential of cosine type and therefore Im Tj = always satisfy the extremal conditions. 
The Im Tj = solutions also decouple the real sector from the imaginary sector completely in 
the mass matrix [7]. Thus we can safely neglect the imaginary components. 



The potential is expanded with respect to (4.2), so 



Defining the following parameters: 

X - = a-t- 6- = ^ C = ~ 27W ^ /2 B =^ (44) 
7lQl 3/2 64V27i^i A i 

the potential takes the approximate form: 

AiWoal ( 2C x ie ~ Xl 



V 



(4.5) 



3 / 2 S.rJ'l 2 \2 



j=2 \ X 1 ^j=2 "j^j 



Here, we have a symmetry of exchanging moduli fields Xj together with the assigned parameters 



.Bj, coming from the symmetry for Tj (i 1) which the full potential (4.1) contains. 

Our purpose here is twofold : to study (1) the effect of increasing the number h ' of Kahler 
moduli and (2) the effect of different distributions for Wq and Ai, on both P(A) and (A). For 
the first case, we shall simply treat Wo, Ai as random parameters with uniform distributions. 
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Figure 10: The probability distribution P(A) of A assuming uniformly distributed Wo,Ai. (A) as a function 
of the number N K of Kahler moduli; the curve is given by (A) - O.OOllOA^ ' 226 ^ - 0268 ^ . 



4.1 Uniformly distributed Wo, A{ 



Now we find the meta-stable minimum of V (4.5) to obtain A and then estimate its probability 



distribution P(A) and (A). Here again we set aj = % = £ = 1 to see just the effect of randomness 
of Wq, Ai on -P(A) and (A). We shall start with uniform distributions for the parameters Wo, A^. 
Here we compare several different distributions with different setups to see how P(A) and (A) 
behave. 

First we assume uniformly distributed random parameters Wo, Ai with range: — 15 < Wo < 
0, < Ai < 1 for convenience. The choice of sign here is for increasing chance of dS solutions. 
Even if we use the different choices, the resultant behavior is unchanged. The calculation cost 
for solving for ti and checking the stability can be hugely reduced if we limit ourselves in the 
parameter region which is valid for meta-stable vacua at A > 0. We follow the analysis of the 



parameter region up to three-moduli in appendix [Dj and further apply the similar limitation 
for the models with more Kahler moduli. The limitation of our interest is the one we obtain 
in the case with three Kahler moduli, but we can expect the similar validity region for B ri — Bj 
plane (i ^ j) as that of B2 — -B3 plane, owing to the symmetry. This limitation hugely reduces 



our calculation costs so we can numerically get the statistical distributions of the potential (4.1) 
even at larger numbers of Kahler moduli. 



In Figure 10, we show plots of the statistical analyses. If we assume uniformly distributed 
Wq, Ai, then the probability distribution -P(A) of A become more sharply peaked as the number 
Nk of Kahler moduli increases. The plot in RHS suggests that the expectation value (A) of A 
is suppressed clS cL function of Nk ■ 

(A) ~ 0.00110iV/- 226 e- a0268 ^, (4.6) 

where we neglected the point at N K = 1 in estimation, since the single modulus model does not 
have off-diagonal components in the Hessian, and therefore we expect that the value at A# = 1 
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Figure 11: The expectation value of the lowest eigenvalue of the physical mass (squared) matrix and the ratio 
of the expectation values. The estimated curve on the LHS goes as (m^ iri ) ~ O.OSlTiV^ 1 - 21 ^ ' 0422 ^ . RHS 
gives the ratio (A) / (w^A as a function of Nk- 

is not necessary to obey the estimated curve from the data among Nk = 2 — 6. 

Although we see the suppression of (A) as increasing Nk, we may worry that the similar 
suppression happens even for mass terms, because the mass matrix is also given as a complicated 
function of the parameters and moduli. We now estimate the physical mass matrix, which is 
the mass matrix for canonically defined fields. The physical mass matrix is related to the 
second derivative of potential d^d^V through the diagonalization matrix of Kahler metric. The 
diagonalization matrix of the Kahler metric is defined to satisfy 

D ki K i3 D% = X k S kl (4.7) 

where Aj are the eigenvalues of Kahler metric and same indeces in LHS are summed over. Then 
the canonical coordinate Yi is defined to be 

dY t = ^i{D T )^dt k . (4.8) 

Now we know the conversion matrix to physical mass matrix from the second derivative matrix: 

2 dt k dti . 
^=d^WU- (4-9) 

Performing the conversion at each stable point, we can calculate the physical mass matrix 
accordingly. 

In the distribution of lowest eigenvalue of mass matrix, we see the more sharply peaked 
behavior as increase of Nk as well as the distribution of A. But as in Figure [TTJ the suppression 
effect is actually stronger than that for (A). The expectation value of the lowest eigenvalue m min 
is estimated as a function of Nk, again neglecting the point at Nk = 1, by 



(m min > ~ 0.0317A^ L21 e- 0422 ^. (4.10) 
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Now we see that the suppression of (m^ in ) is faster than the suppression of (A) because of the 
difference in the exponent. In addition to the similar suppression we have in A, the non-trivial 
off-diagonal components cause the eigenvalue repulsion in the mass matrix when the rank of the 
matrix gets higher. As a result of this, the minimal eigenvalue of mass matrix is further pushed 



down in general. This tendency can be read clearly from the RHS plot in Figure 11 



Once the lowest eigenvalue of mass matrix reaches smaller values, we may worry about the 
cosmological moduli problem 22 24 , implying m min > O(10) TeV ~ 10 _15 Mp[^J In the current 



situation with uniformly distributed Wo, the lowest eigenvalue of mass matrix can be smaller 
than the A. Therefore we encounter the cosmological moduli bound (m^) ~ 10~ 30 where 
Nk ~ 1-35 x 10 3 , earlier than the value for the current cosmological constant (A) ~ 10 -122 
where N K ~ 1.01 x 10 4 . 



4.2 Sharply peaked P{Wq) and uniform distribution P(A') 

Motivated by the analysis in previous subsection, it is quite possible that the parameters W , A,- t 
themselves are given as products or non-trivial functions of random parameters as a result of 
the stabilization of the complex structure and dilaton moduli. As a toy model, let us consider 
a peaked distribution P(Wq) for Wq while Ai remain to have a uniform distribution as before. 



Motivated by the product distribution case (2.2), we simply assume 

1 / 15 \ n ^ 

where again the choice of sign does not affect the result. The number "15" is just a choice 



to increase the chance of getting dS solutions in the numerical simulation. As we see in (2.4), 
the expectation value of Wo is exponentially suppressed in the presence of multiple random 
parameters. 

Here we also assume that the values of s, Ui are given as roughly of order one, and therefore 
negligible. Even though we observe that the s, Ui behave differently in the simplest model 



for complex structure moduli stabilization in section |3.4[ it is possible that the situation will 
be totally different in a different class of models for complex structure moduli stabilization. 
But since we do not have enough information of the other models, here we simply ignore the 
direct effect of s, Ui so that we can simply extract the effect of the combination between peaked 
distribution of Wo and Kahler moduli dynamics. 

In this setup, we solve the model numerically satisfying the stability constraint at positive 



A, up to n = 5 at each Nk = 1, 2, 3. The resultant plots are showed in Figure 12 The estimated 



lr The stringent bound for the reheating temperature T r > 0.7 MeV comes from successful neutrino thermal- 



ization and BBN 28 29 . 
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Figure 12: The expectation values where the distribution P(Wq) has a (n— 1) th power of logarithmic divergence 
at Wo — 0, at Nk = 1 (red, line, circle), Nx — 2 (blue, dashed, square), and Nk = 3 (green, dotted, diamond). 
Each expectation value is suppressed exponentially, though the exponential suppression at Nk = 1 is most 
significant in (A). Note that the first two plots are given as log-plots. 



curves for (A) are 



(A)^ =1 =0.00201n a644 e- a597n 



N K =1 , 

(A) Nk=2 =0.00164n a44 V a616r \ (4.12) 

(A) Nk=3 =0.00127n a409 e-°- 491n . 

When we compare the coefficients of exponent, which is important for larger n, we see that the 
effect of the suppression is most significant around Nx = 2 though the increase of n seems to 
work to make the cosmological constant smaller generically. 

On the other hand, the expectation value of minimal eigenvalue of the mass matrix are 
estimated by 

(mLn) NK ^ =0.0755n°- 72 V - 620 ", 

(mL*)^ =0.0224n°- 422 e-°- 59 -, (4.13) 
«*) NK=3 =0.0115n°- 364 e-°- 448 ". 
We now see that the coefficients in exponent is quite similar to that of the one for (A). Since 



we have the Wq dependence in the coefficient of the potential (4.5) and this coefficient is the 
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common part between potential and mass matrix, the similar suppression is expected even in 



the eigenvalue of mass matrix (see also the last plot in Figure 12). 

Together with the estimated functions above, now we discuss about the cosmological moduli 
problem. When we estimate the values A ~ 10 -122 , Tn^ aia ~ 10~ 30 with the fitted functions 
above, we get the numbers: 



N K 


1 


2 


3 


(A) ~ lO" 122 


n ~ 467 


n ~ 450 


n ~ 564 


«i„> ~ io- 30 


n ~ 113 


n ~ 113 


n rsj 148. 



(4.14) 



The estimated n for the lower bound of moduli mass is always smaller than that for A. Therefore 
we still encounter the cosmological moduli problem. 



4.3 Sharply peaked P(W ) and P{A l ) 

In previous subsection, we considered sharply peaked P(Wq) but uniform distribution for Ai. 
Here we consider the situation where both P(Wq) and P(Ai) are sharply peaked. In this setup, 
the probability distribution functions of the parameters are given by 



P(W ) 
P{Ai) 



15(n- 1)! 
1 



(n 



An- 


15 


w 








T 







n-l 



, -15 < W < 0, 
< Ai < 1. 



(4.15) 



In the numerical calculation, we can get a result even at n = 15 in this setup. This is probably 
related to the preference of small values of Bi = Ai/A\ to stabilize the model at positive A. Now 
since the distribution of Ai prefers smaller values, the resultant E>i are also likely smaller than 
that in case of uniform distribution. 



The plots for the expectation values are given in Figure 13 The estimated functions for (A) 
here are given by 



(A)^ =1 

( A ) Nk =2 
( A ) Nk =3 



=0.00468n 0803 e- L40n , 

^o.oosron - 97 ^- 1 - 49 ", 



,1.48 -1.55n 



=0.00340n 1 

It is interesting that each coefficient of the exponent is below —2 In 2 



(4.16) 



-1.39. Since the 



coefficient appearing in (4.5|) is given as a product of At Wo, we may naively expect an exponential 
suppression of the form e" 



-2raln2 
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Figure 13: The expectation values with peaked distributions P(Wq) and P(Ai) for the parameters Wo,Ai, 
with (n — 1) th power of logarithmic divergence at Wo — and Ai — respectively. The expectation values are 
given for Nk — 1 (red, line, circle), Nk — 2 (blue, dashed, square), and Nk = 3 (green, dotted, diamond). Each 
expectation value is suppressed exponentially, though the exponential suppression at Nk = 1 is most significant 
in (A). Notice that the first two plots are given as log-plots. 

The estimated functions for the lowest eigenvalue of mass matrix become 

Kin)^ =0.175n°- 139 e- 1 - 40 ", 

« in }^ 2 =0.0609n°- 732 e^ (4.17) 



Again, each coefficient of the exponent is below —2 In 2 ~ —1.39. The non-trivial complicated 
moduli couplings push the expectation value down further in addition to the suppression by the 



coefficient. But here, as in the second plot of Figure 13, the expectation value at larger Nk with 
fixed n always stays at lower value, meaning that the off-diagonal components in the Hessian 
tends to push the lowest eigenvalue further down as the rank Nk of the Hessian increases. This 
effect tends to push up the ratio (A) / (m^ in ). 
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Now we are ready to consider the cosmological moduli problem. 



N K 


1 


2 


3 


(A) ~ 1(T 122 


n ~ 197 


n ~ 188 


n ~ 182 


(m 2 ■ ) ~ 10~ 30 

\ mm/ 


n ~ 48 


n ~ 44 


n ~ 42. 



(4.18) 



We see that the required number of n for the lower bound of mass is always smaller than that 
for the cosmological constant A; therefore we cannot avoid the cosmological moduli problem. 



If we compare the result here with the one in previous subsection (4.14), we see that the 
required n here is much smaller. This is nothing but a benefit of the additional peakiness in 
Ai, which is considered to come from the effect of multiple moduli stabilization in the complex 
structure sector. Although so far we have considered up to Nk = 3 in this subsection just for 
simplicity, it is quite possible that we can explain at least the small value of A naturally in terms 
of the variety of string landscape. But we also need to handle the cosmological moduli problem 
simultaneously in some way, which would be quite ubiquitous when we consider multiple moduli 
space because the eigenvalues repulsion gets more serious in the higher rank mass matrix. 



5 Discussions and Remarks 

When there is no preference for a small cosmological constant A, extreme fine-tuning is necessary 
to end up with an exponentially small A, when the Planck mass Mp is the only scale in the 
model. Here we show that simple properties of probability theory can lead to a vanishingly 
small A in some stringy scenarios. In particular, the peaking of the probability distribution 
.P(A) at A = is quite robust. Admittedly, the particular model we study is chosen more for 
its solvability than for its phenomenological properties. Nevertheless, it is a model that have 
been extensively studied in the literature and is non-trivial enough to present some interesting 
features to serve our specific purpose. This possible stringy mechanism for the preference of an 
exponentially small A should be further explored in more realistic stringy scenarios. 

In this paper, we consider the simplest model for the complex structure moduli stabilization, 
in which we only have linear terms of the complex structure moduli U{ in the superpotential 
W Q . We can expect more peaked distributions in w if there are cross couplings among Ui 
in the superpotential. Actually, such cross couplings are quite ubiquitous even in toroidal 
orientifolds |15|. Therefore it would be important to investigate the probability distributions 
P(wo) and -P(A) in the more realistic models to see whether and under what conditions A 
will be naturally vanishingly small. We are encouraged by a preliminary investigation where 
the peaking behavior of P(wq) seems to strengthen when quadratic terms of complex structure 
moduli are introduced into the superpotential. Obviously a more detailed study is required to 
get a full picture of the distribution of A. 
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In this paper, we assume all parameters are physical parameters that include radiative cor- 
rections. In this sense, we are calculating P(A) for the physical A. We are interested in the 
situation that the corrections are small compared to the original inputs and therefore can be 
effectively absorbed in the parameters assigned in the model, which we randomize. If we let 
all the parameters in the model be bare parameters, then the calculated P(A) is obviously for 
the bare A. In general, when the parameters are randomized, the ranges and distributions of 
the bare parameters and that of the physical parameters will be different. Since we have only a 
vague notion of their ranges and distributions, it is important to know when their ranges and 
distributions will make a difference in the resulting physics. As we believe, the peaking behavior 
of P(wq) and P(A) are quite generic for reasonable ranges and distributions, while (|A|) is more 
sensitive to the details. 

Even if the radiative and other corrections contribute significantly to the potential, we expect 
that the peaking behavior is still there since the correction terms will bring in non-trivial moduli 
dependence and the resulting correlation most likely strengthens the peaking behavior in general. 
It is interesting to see if how the peaking behavior is affected by the presence of the other types of 
a' and string loop corrections. Also the backreaction of the a'-correction and non-perturbative 
terms generating the potential for Kahler moduli would be important for complex moduli sector. 



Recently it is argued that the D7 tadpole cancellation condition 30 of the system and the 
holomorphicity of D7-branes in flux compactification restricts the maximal rank of gauge group 
on D7-branes 31 , 32 , which is a candidate to derive non-perturbative terms as a result of the 
gaugino condensation. The explicit solution for dS space was obtained with CP X 1 1 69 even 
in the presence of the constraint for Cj (321. Although we have analyzed statistical behavior 
only by randomizing Wq, At with the other parameters fixed just for simplicity in this paper, 
it would be interesting to see the effect by randomizing the parameter cij within the known 
restrictions/constraints. 

We discuss the mass of the lightest modulus. Its very small mass may pose the cosmological 
moduli problem. This problem may be ameliorated in some multi-field inflationary senarios 
[33 35 based on [36]. We also know that radiative correction to the cosmological mass squared 
goes like M| USY where Msusy is the SUSY breaking scale, so the cosmological moduli problem 
may not be as serious as one expects. In any case, a better understanding of this problem is 
important if we like to take this stringy mechanism for a vanishingly small A seriously. If we 
impose the present cosmological moduli constraint, we can bring A down only a limited number 
of orders of magnitude. It is entirely possible that the stringy mechanism discussed in this paper 
provides only a partial solution to the A problem. 
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A Some toy models 

To get some feeling on the peaking property for the probability distribution -P(A) of A, consider 
the single field polynomial model. At least three terms with different powers are required to 



achieve the metastable vacua with a positive A 37 40 . Since we are interested in only the 



meta-stable solution, we may consider a model of a simple minimal form: 

V(0) = a<j> - \f + ^ (A.l) 

where all parameters are positively defined for simplicity, and have the scale Mp. Now, we 

impose the stability d|l/| mm > at the extremal points given by d^V = 0. These conditions 
imply 

A = Vb 2 - 2ac > 0, <p min = . (A.2) 

c 

Plugging back to the potential, we get 

A ^ VU. = (t + A) 6 2 c ( r 2A) . (A.3) 



So we see that the sum form (A.l) is converted to the product form (A.3). 

If we impose the positivity of the potential A > 0, then the parameter region is further 
restricted by 

where the lower bound is given by the stability while the upper bound is for the positivity. 
Whenever we satisfy the region of the combined parameter above, there exists a minimum with 
positive A. Since we are interested in the region with small As, let us further focus on the region 
around C p ~ 4/3. Then the potential becomes 



. 9ab /4 » 

A~--(--C p ). (A.5) 



Next we analyze the distribution of the approximate minima (A.5). To make the argument 



simple, we here fix c = 1 and take a, b as random parameters obeying uniform distribution 
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P(ln(A)) 

0.4 - 




ln(A) 



Figure 14: The comparison of the analytical (the curve) and the numerical (the histogram) distributions P(A) 



vs A (LHS) and P(lnA) vs In A (RHS) in the polynomial model (A.l) 



within a range < a,b < 1. Then the probability distribution function is estimated to be 



P(A) = Ng 1 J dC p J dadb5 



(A.6) 



where iVo is the normalization constant. Performing the last two integrals for the two (^-functions, 
being careful with the integration regions, we get 

1 „ 12 



P(A) 



No 1 



dC n 



3C p (4/3 - C p 



1<C„< 

p ~ 9 + 16A 



(A.7) 



The constraint for C p also suggests < A < 3/16. The normalization constant Nq is calculated 
such that / dAP(A) = 1, which gives A" = 1/24, so 



P(A) 



16 



In 



16A 



(A.8) 



where P(A) is divergent at A = 0. This is nothing but the outcome of the product form of the 



coefficient obtained in (A.5). If we allow the parameter c to have a distribution (even including 
zero), the divergence at A = is still present (see the argument for ratio distributions, e.g. 
in However, P(A) will have a longer tail at large A. 



Although we have analyzed the approximate potential (A.5 ), we can easily check numerically 



the divergent behavior (A.8) at the minimum (A. 3) in the full form of the potential. In Figure 



14, the analytical function (A.8) perfectly agrees with the numerical histogram especially at 



lower values, given uniformly distributed random values for < a, b < 1 while fixing c = 1. 
Next, consider a simple multi-field case, 

v((f)i) = Ai4>i + Bij4>i4>j 

We see that 



V m 



-A T B- 1 A 
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which is a sum of terms for the multi-field cases. Clearly the Py coupling alone is not enough 
for l^nin to take the product form. Higher couplings for tighter interactions among the fields are 
necessary to bring V min to the product form. 

Although we do not have an analytic form for V m i n for more complicated V(<f), it is easy to 
convince oneself (and check numerically) that a divergent P(A) at A = for meta-stable V m i n 
is quite typical for V in which all terms are tightly coupled to each other. 

For more complicated V ((/)), V m i n takes a more complicated form. We are looking for 2 
features : (1) whether the probability distribution P(A) has a peaking (maybe even divergent) 
behavior at A = 0, and the peakiness becomes more pronounced as the number of fields and their 
couplings (parameters treated as random variables) increase; and (2) whether the expectation 
value of the magnitude (A) is dropping exponentially fast as the number of fields increases. 



B -P(A) in the single Kahler modulus case 

We show here the way to estimate the probability distribution P(A) in the single Kahler modulus 
case discussed in section 12.21 Let us define 

Wl = -W , w 2 = A u c = n ^= 64 ^C. (B.l) 

w 2 27 

and set the value a\ — 71 = £ = 1. Neglecting the overall coefficient, the model is simplified to 
be 

z = wiw 2 {c - c ), c < c = — < ci, (B.2) 

U> 2 

where Co ~ 12.2 and c\ ~ 13.0. The cosmological constant is related with this z through the 
relation A = 3^/(2500^). 

Now the probability distribution of z is determined by 

P( z ) = Nq 1 J 1 dc dw 1 dw 2 5{w 1 w 2 {c-c )) 6 - cj , (B.3) 

where A^o is the normalization constant. Since this is the constrained integrations by two delta 
functions, the integrated region is further constrained. The last two integrations can be achieved 
easily and we get 

p(z)= j N o l !dc^ for0<z<c(c-c ), c < 1, 
" S N o l !dc^ for0<z<^, c>l. 
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Since the region of c of interest is 1 < Co < c < ci, we can integrate over c and 



P(z) =N 1 



N 



Cl 



dc 



l-z 



2c(c - c ) ci - c 



In 



Ci - c 



Ci2 



<fc / dc — r 

2c(c — c J 



Ci - c 



(B.5) 



The probability distribution of A can be estimated via the relation P(A) dA = P{z) dz by. 

3 ci - c 



2500V5 g, ln 



Ci — Cq 



2500^ ciA 



(B.6) 



which is (2.12). 



C Special symmetric case for s and kjq 

Some of the manifolds have special symmetries so it is interesting to see the impact of them on 
the value of s and the probability distribution P(wq). 

Let Ti = di/bi and r = C2/C1. Let us consider the special situation, 

(Cl) 



r = ri = T2 



while r 7^ tq. Following from (|3.7|), all complex structure moduli Ui are equal. Then the solution 

(C2) 



of (3.10) is given by 



1 



2rr 



(h 2 ' 1 - l)(r - r ) ± ^(/i 2 ' 1 - l) 2 (r-r ) 2 + 4rr 



which, for large h * yields, for r > r , 

s ~ (/i 2 ' 1 - 1) 



r — To 



1 



or 



(C.3) 



rro (/i 2 ' 1 — l)(r — ro) 

in which case, r must be close to ro and tq > r for s to be large and positive. In this case, we 
now have 

Cl 



W, 



I min 



h 2 



. ( 1 — sr \ 



(C.4) 



It is easy to see that the strong peaking behavior of P(wq) is absent here. 



At the symmetric point, &j = di or r = = 1, ( |C.3[ ) reduces to 



(h 2 ' 1 - 1) 



Cl 

c 2 



(C.5) 



Recall that we like s to be large and positive, say s ~ 10, so we need ci > c 2 . 

The following picture emerges. As we move towards the symmetric points, the probability 
distribution P(wq) loses its peaking behavior at A = 0. This is simply because (the random) 
parameter space collapses at the symmetric points. 



36 




Figure 15: The validity region in B2 — C plane for meta-stable dS vacua with two moduli. The x-axis is B2 
and the y-axis is C. The contours inside the region describe level lines of U, not the magnitude of A. 

D Stabilization in two and three Kahler moduli models 

For the analysis of the potential for multi-Kahler moduli in section |4j we elaborate on the 
stability analysis in the cases of two and three Kahler moduli. Here we just focus on the bracket 



part of (4.5) since the coefficient part does not touch the stability. For convenience, we further 
define 

_ 2C x 1 e~ Xl ^ BiXie~ Xi 

= 9(*? /2 - E i=2 ^f )* " W 2 - ~ h W 2 - e~^?¥' (m) 

D.l Two Kahler moduli 

Let us start with the two moduli case. In the region (S 2 < 1, the stabilization of x\ is little 
changed from that of the single Kahler case. Therefore the moduli stabilization for X2 requires 



small B2 so the third term in (4.5) is comparable with the other terms proportional to 82 ■ The 
situation begins to differ from the single Kahler case as 62 increases. But since the volume is 
defined to be positive and we expect the fixed moduli value of Xi to be comparable with x%, 62 
should not be too large. So hereafter we take 8 2 — 1 just for simplicity. For stability, we need 



to be careful about the positivity of mass matrix. According to Sylvester criteria (see e.g. 41 
which is also used in 42 ), we consider positivity of d^U and the determinant of the 2x2 matrix 



d Xi d x .U , simultaneously with the positivity of the extremal points. 
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The extremal conditions d Xi U = are solved by 
C - 



xl(x 2 - 1) 

x (2x1 + x\){x 2 - 1) - xl /2 x 3 2 /2 (x2 - 4) + (xxxl - x\)(x 2 + 2) - x\ ,2 x^ 2 (2x 2 + 1) 

(D.2) 

If we take x 2 — >■ 0, the conditions recover that for the single Kahler case. Together with the 



positive mass squared conditions above, Figure 15 shows the allowed parameters region with 



respect to C, B 2 for meta-stable dS vacua. In the system above, the positivity constraint of the 



determinant controls the upper bound of the parameter region in Figure 15, while the lower 
bound is given by the positivity constraint U > 0. The contours inside the region describe level 
lines of U, not the magnitude of A. On the other hand, A = at the top boundary curve, and 
A increases as we move down away from this boundary. The numerical values of the parameters 
in the region are roughly given by 

3-650 < C < 4.498, < B 2 < 0.1769. (D.3) 
For instance, the lower bound of the parameters requires the potential minimum at 

Xl = 2.5, x 2 = 0, (D.4) 
while the upper bound of the parameters suggests the minimum at around 

xi ~ 3.881, x 2 ~ 0.7085. (D.5) 

Note that the value of U increases gradually with C, similar to the single Kahler modulus case. 

Note that the validity range of the C parameter is different from the one obtained for single 
Kahler moduli. The result for single Kahler modulus case is (7j 



3.650 < C < 3.887, B 2 = 0. (D.6) 

Therefore we see that the additional modulus dependence alleviates the constraint for C param- 
eter, but still restricts B 2 small. We may understand the smallness of B 2 as a consequence of the 



positivity of the volume V defined in (4.1). The required smallness of ^ i=2 liU/liti affects the 
combinations AijA\ to stabilize the potential in entire Kahler moduli space. Therefore A{jA\ 
cannot be large and are about the same order of Xijx\ (or the combination of 5iXi/x) at the dS 
minima. 
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Figure 16: The parameter region in the B 2 — C — B 3 plot shows the stable minima at positive A with three 
moduli. The last two figures display the data on a projected plane. 



D.2 Three Kahler moduli 

Now let us proceed to the model with three Kahler moduli case. Since it is difficult to draw figures 
with continuous surfaces, here we use the numerical data given random values for parameters 
assigned to the model. Similar to the two moduli case, we set 5 2 = S 3 = 1 just for simplicity. 

Given the values for parameters C,B 2 ,B 3 obeying uniform distribution, we get the region 



of parameter sets which satisfy the stability condition at positive A, shown in Figure 16 All 
points here satisfy the stability constraint at extremal points with A > 0. We see that the 
stable parameter region for three moduli case includes the stable region for two moduli case as 
a boundary. The restricted ranges for the parameters are given by 

3-650 < C < 4.498, < B 2 < 0.1769, < B 3 < 0.1769. (D.7) 

The parameters are correlated so they are further constrained. 



D.3 Checking the validity of the approximate potential 



So far we have seen what happens in the approximate potential (D.l). Here let us compare the 



resultant parameter space with that of the full potential (4.1). The full potential is no longer 



written down in terms of the combined parameters (4.4). Therefore we give values for Wo,Ai 
and set a>i = 7, = £ = 1 for simplicity, and solve for moduli fields ti. 

The full potential for two Kahler moduli suggests the stable parameter region for C, B 2 as 



in LHS of Figure 17, where 



3-95 < C < 4.87, < B 2 < 0.193. 



(D.8) 



Although the boundary values for parameter region are changed a little from the approximate 
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Figure 17: The stable parameter region for dS vacua in the full potential (4.1) for the two and three moduli 
case. 



ones (D.3) for the approximate potential (4.5), the essential feature is unchanged, including the 



shape of the region which is crucial for the distribution of A. 

In the three moduli model, the full potential also suggests the similar parameter region for 
stability to that for the approximate potential. The stability condition at positive A is satisfied 



if the parameters are within the region in RHS of Figure 17, suggesting roughly 



3.95<C<4.87, 0<£? 2 < 0.193, < B 3 < 0.193. 



(D.9) 



Again, though the values are changed a little from (D.7), the shape of stable region is the same 
as before. 

From the examples of the analysis of the full potential above, we see the basic agreements 
between the approximate potential and the full potential, at least for the distribution analysis 
of small A. This approximation is quite helpful if we like to go beyond three Kahler moduli. 



Thus, as in single modulus case in section 2.2, we may use the approximate potential (4.5) to 



understand the statistical properties of the potential (4.1) for the general multi-Kahler moduli 
cases studied in section HI 
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